Entanglement, fidelity and topological entropy in a quantum phase transition to topological order 
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We present a numerical study of a quantum phase transition from a spin-polarized to a topologically ordered 
phase in a system of spin- 1/2 particles on a torus. We demonstrate that this non-symmetry-breaking topological 
quantum phase transition (TOQPT) is of second order. The transition is analyzed via the ground state energy and 
fidelity, block entanglement, Wilson loops, and the recently proposed topological entropy. Only the topological 
entropy distinguishes the TOQPT from a standard QPT, and remarkably, does so already for small system sizes. 
Thus the topological entropy serves as a proper order parameter. We demonstrate that our conclusions are robust 
under the addition of random perturbations, not only in the topological phase, but also in the spin polarized phase 
and even at the critical point. 
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I. INTRODUCTION 

A quantum phase transition (QPT) occurs when the or- 
der parameter of a quantum system becomes discontinuous 
or singular—. This is associated with a drastic change of the 
ground state wave function. Unlike classical phase transi- 
tions, QPTs occur at T = and thus are not driven by ther- 
mal fluctuations. Instead, quantum fluctuations are capable of 
changing the internal order of a system and cause the tran- 
sition. When a quantum Hamiltonian H{X), which depends 
smoothly on external parameters A, approaches a quantum 
critical point Ac from a gapped phase, the gap A above the 
ground state closes, and the critical system has gapless exci- 
tations. This corresponds to a continuous, second order QPT. 

Here, we consider a QPT from a spin polarized to a topo- 
logically ordered phase: a topological quantum phase tran- 
sition (TOQPT). The internal order that characterizes topo- 
logically ordered phases cannot be explained by the standard 
Ginzburg-Landau theory of symmetry breaking and local or- 
der parameters. Instead, it requires the notion of Topological 
Order (TO)^. TO manifests itself in a ground state degener- 
acy which depends on the topology of the physical system, 
and it is robust against arbitrary local perturbations-'. This ro- 
bustness is at the root of topological quantum computation, 
i.e., the ground state degeneracy can be used as a robust mem- 
ory, and the topological interactions among the quasi-particles 
can be used to construct robust logic gates'*-^. On the other 
hand, to what extent a TOQPT is affected by perturbations is 
a problem that has only very recently been addressed''-^, and 
is a focus of this work. Moreover, the classification of TO is 
still an open question. Ground state degeneracy, quasiparti- 
cle statistics and edge states, all measure and detect TO but 
do not suffice to give a full description. Tools from quan- 
tum information theory, specifically entanglement^-^ and the 
ground state fidelity'", have recently been widely exploited to 
characterize QPTs. To date, all the QPTs studied with these 
tools have been of the usual symmetry breaking type. Here 
we apply them to the transition from a spin-polarized phase to 
a TO phase, and find that they are universal in the sense that 
they detect this transition. However, these tools do not suf- 



fice to distinguish a symmetry breaking QPT from a TOQPT. 
Recently, the new concept of "topological entropy" 5top was 
introduced". The topological entropy vanishes in the ther- 
modynamic limit for a normal state, whereas S'top ^ for 
a TO state. Therefore, S'top can serve as an order parame- 
ter Moreover, TO is not only a property of infinite systems, 
and an important question that was left open in Refs." is the 
behavior of S'top for finite systems. Here we shed light on 
this question by presenting finite-system calculations of Stop. 
We report that Stop changes abruptly at the critical point of 
a phase transition between phases with and without TO, even 
for very small systems. It is thus an excellent discriminator 
between the absence and presence of TO, and moreover. Stop 
is capable of detecting a TOQPT. 



Specifically, we present an exact time-dependent numerical 
study of a TOQPT, introduced in Ref.^, from a spin-polarized 
phase to a TO phase, for both the ideal model and the model in 
presence of an external perturbation. Our results are the fol- 
lowing: (i) standard QPT detectors (derivative of the ground 
state energy', entanglement of a subsystem with the remain- 
der of the latticd^iS, ground state fidelitjiiS), are all singular at 
the critical point of the TOQPT, thus confirming that this is 
indeed a QPT. Ground state fidelity and block entanglement 
are thus capable of dealing also with non symmetry breaking 
QPTs. (ii) Stop detects the TOQPT in a very sharp manner 
already for small system sizes. It also detects TO better than 
other non-local order parameters, in particular the expecta- 
tion value of Wilson loops. It is therefore appropriate for the 
detection and characterization of TOQPTs and for studying 
TO. These results complement and strengthen the conclusions 
of Ref. ' ' . (iii) Adiabatic evolution can initialize topological 
quantum memory faithfully: even in the presence of pertur- 
bations the coupling to other topological sectors and excited 
states is negligible, (iv) This robustness extends to the entire 
topological phase, and even to the critical point itself. Pertur- 
bations do not affect the nature of the TOQPT either. 
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FIG. 1: (Color online) A square lattice with 32 spins. The spin de- 
grees of freedom are placed on the vertices. The red dashed lines 
tf , t2 are the incontractible loops around the torus. The product 
B7 ■ B8 ■ Bll denotes the loop operator drawn in red. All the spins 
on the vertices crossed by a loop are flipped. The region AuBuC is 
a ring containing eight spins, used in computing Stop- For the lattice 
of 32 spins, the ring has diameter R = 2 and width r — 1. 



II. PRELIMINARIES 

Consider a square lattice L with periodic boundary con- 
ditions (torus) and with n spin- 1/2 degrees of freedom oc- 
cupying its vertices. The Hilbert space is given by 7i = 
span{|0), ll)}**", where |0) and |1) are the ± eigenvectors 
of the Pauli cr^ matrix. As shown in Fig. [l] the n pla- 
quettes can be partitioned into two sub-lattices, denoted by 
different colors. Following Kitaev**, we associate with ev- 
ery white plaquette p an operator Bp = Yij^dp '^j ^^P^ 
all spins along the boundary of p. A "closed string opera- 
tor" is a product of plaquette operators Bp that flips all spins 
around a loop (or around a loop net). The "group of closed 
strings" X is the group of products of plaquettes Bp. Sim- 
ilarly, with every pink plaquette s, we associate an operator 
As = TTocc which counts if there is an even or odd num- 
ber of flipped spins around the plaquette s. Kitaev's toric 
code Hamiltonian'* is then given by Hjjg — —Uj^s ~ 
g'}2,pBp = Hjj + Hg, which realizes a Z2 lattice gauge 
theory in the limit U ^ 00. The ground state is an equal 
superposition of all closed strings (loops) acting on the spin 
polarized state |vac) = |0)i ® ... ® |0)„ - it is in a string- 
condensed phase. The ground state manifold is given by 
£ = span{|X|-5(if)'(tfp^^^^x|vac); e {0,1}}, 
which is fourfold degenerate'^. The if 2^ Aip the spins 
along an incontractible loop around the torus (See FiglT), tak- 
ing a vector in C to an orthogonal one in the same manifold 
because they commute with Hu g. On a lattice on a Riemann 
surface of genus g, there are 2g incontractible loops {t^ }^li, 
and therefore C is 2^3-fold degenerateiii (for a torus g = 1). 

The Model and the QPT. — Now consider the following 
time-dependent Hamiltonian, introduced in'' as a model for 



a TOQPT: 

Ho{r)^Hu + rHg + {l-T)H^, (1) 

where = X^r^i o-^, t = t/T € [0, 1], and T is the total 
time. The non-degenerate ground state of H{Q) — Hu + i?^ 
is the spin polarized state |vac) which is the vacuum of the 
strings. The term (1 — t)H^ acts as a tension for the strings, 
whereas rHg causes the strings to fluctuate. As t increases, 
the string fluctuations increase while the loop tension de- 
creases. For a critical value of A = Tg/{1 — r)^, and in the 
thermodynamic limit, a continuous QPT occurs to a TO phase 
of string condensation. This QPT is not symmetry breaking, 
i.e., is a TOQPT. As argued in Ref.'', provided T > 1/A,nin 
(the minimum gap, as a function of t, between the ground 
state and the first excited state) evolution according to H{t) is 
an adiabatic preparation mechanism of a TO state: one of the 
degenerate ground states of Kitaev's toric code model-. 
Ref.^ showed that Amin ~ l/\Ai- H{t) can be mapped 
onto an Ising model in a transverse field, which is known 
to have a second order QPT— (see also^). However, in this 
work we do not resort to such a mapping, because it is non- 
local and does not preserve entanglement measures. Instead, 
we numerically study H{t) for r G [0, 1] in At = .01 in- 
crements on lattices Ln with n = {8, 18, 32} spins, and set 
U = 100,^ = 5 = 1. The computational methods used here 
are (i) the Housholder algorithm.'^ for the full diagonalization 
(all eigenstates) of L8, and (ii) a modified Lanczos method — 
to obtain the low-energy sectors of L18 and L32. We observe 
that for all r G [0,1] the ground state comprises only closed 
strings. Since this is the case for every finite system size, and 
in order to reduce computation cost, we diagonalize L32 only 
in the relevant symmetry subspaces, defined by the constraint 
^.l^)^n,e.^|IV'> = lV'>,Vs. 



III. THE PERTURBED MODEL 

To test the robustness of the TOQPT, we also studied the 
perturbed model given by 

n 

H{t) = Ho{t) + Ho{t) + J2 (/»"(J>J + h'{j)<y]) 

(2) 

The perturbation V is random with (j) and /i^ (j) uniformly 
distributed in [—0.2,0.2] and [— P, P], respectively, with the 
magnitude P variable in our calculations below. We carried 
out calculations for L8 (time-dependent) and L18 (ground 
state only). These were averaged over random realizations of 
V, and included the full Hilbert space as V disrupts the sym- 
metry A s\ip) = IV')- The z-component of the perturbation is 
expected to have a small effect as it only slightly modifies the 
term for r < Tc, while for t > Tc TO dominates and ten- 
sion effects are suppressed. Our calculations confirmed this, 
and hence Figs.|2]|7]show the results for = 0. 
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FIG. 2: (Color online) Fidelity between the time-dependent solution 
of the Schrodinger equation and the adiabatic state, for different val- 
ues of the total evolution time: T = 20, 40, 60. (a) The unperturbed 
model for L18. The evolution is adiabatic for T = 60. Note that 
the drop in adiabaticity is a precursor of the QPT. (b) L8: fidelity 
in both the ideal and perturbed (P=l) cases. The perturbed model is 
indistinguishable from the ideal one. 



IV. ADIABATIC EVOLUTION 

We numerically simulated the time evolution from the fully 
polarized state at t = to the string-condensed phase at 
T = 1. The possibility of preparation of topological order 
via such evolution has been studied theoretically in Ref. ,6. A 
crucial point is to show that the adiabatic time depends on the 
minimum gap that marks the phase transition (and that is poly- 
nomially small in the number of spins), and not on the expo- 
nentially small splitting of the ground state in the topological 
phase. To this end, one must show that transitions to other 
topological sectors are forbidden and protected by topology.* 
The initial wave function is the exactly known ground state 
of H{t = 0). This state is then used as the seed to com- 
pute the ground state of H(Ar). After iteration, this state 
is in turn used as the seed for H{2At), etc. We can esti- 
mate to what extent the evolution is adiabatic by numerically 
solving, for LIS, the time dependent Schrodinger equation 
Hi/j{t) = i-ip{T) for different values of the total evolution 
time T. This is shown in Fig.|2|a), where we plot the fidelity 
between the time evolved wave function V'(t) and the instan- 
taneous ground state: JFad = l(''/'(''")IV'o('''))|- Moreover, we 
compute JFad also for the perturbed model, but the largest lat- 
tice for which we can do this is L8. Fig. 2b shows clearly that 
for P — 1 the perturbation does not change the time-evolved 
state. Significant effects start at P — 2 (not shown). We also 
find that the overlap between the evolved wave function ■i/j{t) 
and the other sectors (ii)*(^2)"' lV'o(''')) is of order 10^'^ for 
every ^ (0, 0) and value of T tested. This is numeri- 
cal evidence for the argument that time evolution will always 
keep the instantaneous eigenstate within a topological sector, 
even in presence of perturbations*. Thus the relevant gap for 
adiabatic evolution is that to the other closed string excited 
states, which implies that the evolution into the TO sector can 
be used to prepare a topological quantum memory*. Hence- 
forth we work only in the sector (i = 0, j = 0), into which the 



system is initialized as the unique ground state of (r = 0). 



V. DETECTING THE QPT WITH STANDARD MEASURES 

To check that the transition from magnetic order to TO is 
indeed a QPT, we first computed the energy per particle of the 
ground state for L8, L18, L32, and its second derivative. As 
seen in Fig. [Ua), the latter develops a singularity as system 
size increases, signaling a second order QPT with a critical 
point atr ~ 0. 71, corresponding to a ratio ^/g ~ 0.41. This is 
in good agreement with the analytical stud}^, which obtained 
(in the thermodynamic limit) g ^ .44, even if this model is 
only asymptotically equivalent to the toric code in a magnetic 
field, in the small field limit. On the other hand, Ref.^ found 
^/g ~ 0.33, using a mapping to the classical 3D Ising model. 
In Fig. we show the block entanglement between four 
spins in a small loop {BW, Fig.[T]) and the rest of the lattice, as 
measured by the von Neumann entropy. In agreement with the 
general theory^, the derivative of the entanglement diverges at 
the critical point for a second order QPT. 

A new interesting alternative characterization of QPTs can 
be given in terms of the scaling in the fidelity Ti\r(j) = 
— Ar))| between two different ground states'". At 
a quantum phase transition, the fidelity should scale to zero 
superextensively. Previous wor k'°'^° has shown that the fi- 
delity criterion is valid for generic symmetry breaking second 
order QPTs. Nevertheless, the fidelity criterion is not strictly 
local, so one would like to know whether it detects the QPT 
to a topologically ordered state. The results are shown in 
Fig. (He). The fidelity drop criterion indeed also detects the 
QPT. Figures [3ja)-(c) also show the result for the perturbed 
model. 

By looking at the behavior of the transition in the presence 
of perturbations, we can safely conclude that the QPT is un- 
affected by the perturbation for P < 10, namely the value of 
Tc and the magnitude of the fidelity drop remain unchanged. 
In Fig. [3jd), we plot the overlap between the perturbed and 
unperturbed ground state. The drop in this quantity also sig- 
nals the QPT, showing that the system is most sensitive to 
perturbations at the critical point (see also Ref.'^). Interest- 
ingly, in contrast to the robustness of the entanglement and 
J^i\t{j), the perturbed and unperturbed ground states differ 
significantly already for P > 2. The results in Fig. [3] thus 
allow us to infer unambiguously that there is indeed a sec- 
ond order QPT in the adiabatic dynamics generated by H{t). 
However, none of the quantities shown in Fig. [3] is explic- 
itly designed to detect topological features, and hence these 
quantities are incapable of distinguishing between a symme- 
try breaking QPT and a TOQPT. 



VI. CHARACTERIZING THE TOPOLOGICAL PHASE 

The spin-polarized regime for t < Tc is characterized by a 
finite magnetization. On the other hand, the topologically or- 
dered phase T > Tc does not admit a local order parameter'-. 
The topologically ordered phase is a string condensed phase 
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FIG. 3: (Color online) QPT detectors for L8, L18, L32, for the un- 
perturbed and perturbed model. All graphs show strong resilience of 
the model and its QPT against perturbations: (a) Second derivative 
of E{t), diverging for Tc ^ .7. The QPT is thus second order, (b) 
Derivative of the von Neumann entropy, measuring the entanglement 
of a plaquette with the rest of the lattice. Its divergence at critical- 
ity also signals a second order QPT. The perturbation has no effect 
for P = 20 (triangles indistinguishable from circles) but is visible 
for P = 40. (c) Ground state fidelity J-{t): the fidelity drop at the 
critical point signals a QPT, associated with a drastic change in the 
properties of the ground state, (d) Overlap between the perturbed 
and the ideal ground state. The clearly visible susceptibility to the 
perturbation at the critical point also signals the QPT. 
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FIG. 4: (Color online) Expectation value of Wilson loop operators of 
increasing size for L32 The expectation value of the loop operators 
starts to increase at Tc, more steeply so for the largest loops, indicat- 
ing that this observable can be used to detect the TOQPT for large 
systems. 
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FIG. 5: (Color online) von Neumann entropy for a plaquette of spins 
and Stop for L32 with an Ising Hamiltonian in a transverse field. 
Note the different vertical scales. For a system without topological 
order, 5top is always ~ 0. (the small bump is a finite-size effect). 



and an effective Z2 local gauge theory and thus the observ- 
ables must be gauge invariant quantities. These quantities 
are the Wilson loops. In this theory, we make a Wilson loop 
^^(2) [-^] of the x{z) type by drawing a closed string 7 on the 
lattice, and operating with (T^(cr^) on all the spins encountered 
by the loop. In the polarized phase, the tension is high and it is 
difficult to create large loops. The expectation value of loops 
decays with the area enclosed by the loop. In the topologically 
ordered phase, large loops are less costly and their expectation 
value only decays at most with the perimeter of the loop. The 
phase transition is of the confinement/deconfinementtype. We 
can write any (contractible) Wilson loop as the product of 
some plaquette operator: M^^(^)[7'] — J^^.^^ i?j,In particu- 
lar at the point r = 1 when the model is the exact toric code, 
the expectation value of Wilson loops is [7] |) = 1 for 

every loop 7, independently of its size. Of course, large loops 
are highly non-local observables. We have computed the ex- 
pectation value of Wilson loop operators of increasing size as 
a function of r. As Fig. |4] shows, the expectation values of 
large loops vanish in the spin-polarized phase, and increase 
exponentially in the TO phase. However, in the limit of infi- 
nite length, Wilson loops are not observables of the pure gauge 
theory^l and cannot be measured. 

Nevertheless, topological order reveals itself in the way the 
ground state is entangled. If we compute the von Neumann 
entropy for a region with perimeter L, the entanglement en- 
tropy will be = L — 1 in the topological phase - see Fig. |6] 
The spin polarized phase is not entangled. We see that there is 
a finite correction of —1 to the boundary law for the entangle- 
ment, which is due to the presence of topological orderi^^ii. 
Therefore we can consider as an alternative non local order 
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FIG. 6: (Color online) S'top for L18 and L32, and von Neumann 
entropy for L32, for the ideal and perturbed model. 5*^ assumes the 
value l — l = 3 in the entire TO phase, where Z — 1 is the exact value 
of Sv for the pure Kitaev model (r = 1) and Z = 4 is the length of 
the border of a plaquette. Stop is zero in the spin-polarized phase and 
quickly reaches unity in the TO phase. Also the topological character 
of the QPT and of the phases is resilient to perturbations. 
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where are the entanglement entropies associated with four 
cuts cr = {AU B U C, AU C,AU B,A}, as depicted in 
Fig. [U We computed 5top(T) in the instantaneous ground 
state \iP{t)) for LIS and L32 (L8 is too small) in the ideal 
model and for L18 in the perturbed model - see Figs. 151617] 
In the spin-polarized phase, even for finite systems, S'top = 
and it becomes different from zero only in the vicinity of the 
critical point, after which it rapidly reaches 1 (as predicted in 
the thermodynamic limit in Ref."). To test whether S'top can 
discriminate between symmetry breaking QPTs and TOQPTs, 
we show in in Fig. |5] the behavior of block entanglement and 
S'top for a quantum Ising model in 2D. This model admits a 
QPT between a paramagnetic and magnetically ordered phase, 
which is symmetry breaking. Block entanglement detects the 
critical point sharply, while Stop does not (note the different 
scales on the left and right vertical axes). The small non-zero 
value of S'top is a finite size effect. 

The block entropy in Fig. |6] shows that the state is already 
rather entangled in the spin-polarized region, whereas S'top is 
almost zero before the transition to TO occurs. Note that the 
block entanglement at the critical point is bounded from above 
by the final-state entanglement (r — 1), which obeys the area 
law. This is an example of the fact that in 2D, critical systems 
do not need to violate the area law as in ID. The useful fea- 
ture of S'top is not only that it can be used in order to locate 
the critical point [Fig.|7l, but also that it allows one to under- 
stand the type of QPT (symmetry breaking or TO). Remark- 
ably, Figs. l5l7l show that S'top has these properties already for 
finite and very small systems. The accuracy of the finite-size 
S'top at the limit points r = 0, 1 is due to the fact that there the 
correlations are exactly zero-ranged. This, however, is not the 
case for intermediate r, especially near the QPT, so how Stop 
works as an order parameter, and how sharply its derivative 
detects the QPT, are rather non-trivial. 

In the presence of the perturbation h^{j), which tends to 
destroy the loop structure, S'top detects the TOQPT up to 
the value P ^ 25, after which a transition occurs: see 
Fig. iTfinset). Overall, Figs. |6|7| show that the robustness of 
TO against perturbations is a feature of the whole topologi- 
cal phase and not only of the analytically solvable model at 
T ^ 1. Finally, we note another remarkable fact: setting the 
x-perturbation V to zero, and moving backward in time from 
T = 1, we can view also the tension term as a perturbation. 
This is due to the fact that the toric code is symmetric under 
the exchange x ^ z in the spin components. The flatness of 
S'top in Fig.|6](squares and circles) shows the robustness of the 
topological phase against this perturbation (see also Refi^). 



FIG. 7: (Color online) First derivative of Stop for L18 and L32, 
for the ideal and perturbed model, diverging at Tc thus signaling the 
TOQPT; (inset) derivative and full width at half maximum of Stop at 
Tc as function of perturbation strength P. Stop remains robust up to 
P ~ 25. 



VII. CONCLUSIONS 

We have presented a comprehensive numerical study of a 
TOQPT. Our results show, using a variety of previously pro- 
posed QPT detectors, that this is a second order transition. 
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Unlike the other detectors, the topological entropy S'top is 
capable of distinguishing this TOPQT from a standard one, 
already for small lattices. Strikingly, the model and its TO- 
QPT are highly robust against random perturbations not only 
deep inside the topological phase, where the gap protects the 
ground state from perturbations, but - even more surprisingly 
- at the gapless critical point. This phenomenon requires fur- 
ther investigation to be properly understood. Moreover, S'top 
detects the TOQPT for perturbations of strength up to 20% of 



the strongest couplings. Of course finite-size effects can be 
important, but it is not possible at present to compute S'top ex- 
actly without direct diagonalization, and this poses limits on 
the maximum size of systems that can be studied. 
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